Gas and stellar dynamics in NGC 1068. Probing the galactic gravitational potential

We present Sauron 2D spectrography of the central 1.5 kpc of the nearby Sey2 galaxy NGC1068, encompassing the well-known NIR inner bar. We have successively disentangled the respective contributions of the ionized gas and stars, thus deriving their 2D distribution and kinematics. The [OIII] and Hbeta emission lines exhibit very different spatial distribution and kinematics, the latter following inner spiral arms with clumps associated with star formation. Strong inwards streaming motions are observed in both the Hbeta and [OIII] kinematics. The stellar kinematics also exhibit clear signatures of a non-axisymmetric tumbling potential, with a twist in both the velocity and h3 fields. We re-examined the long-slit data of Shapiro et al (2003) using pPXF: a strong decoupling of h3 is revealed, and the central decrease in h4 hinted in the Sauron data is confirmed. These data also suggest that NGC1068 is a good candidate for a so-called sigma-drop. We confirm the possible presence of two pattern speeds. We also examine the stellar kinematics of bars formed in N-body+SPH simulations built from axisymmetric initial conditions. These successfully reproduce a number of properties observed in the 2D kinematics of NGC1068, and the long-slit data, showing that the kinematic signature of the NIR bar is imprinted in the stellar kinematics. The remaining differences between the models and the observed properties are mostly due to the exclusion of star formation and the lack of the primary large-scale oval/bar in the simulations. These models suggest that the inner bar could drive a significant amount of gas down to a scale of ~300 pc. This is consistent with the interpretation of the sigma-drop in NGC1068 being the result of central gas accretion followed by an episode of star formation.Comment: accepted for publication in MNRAS, 20 pages, 17 figures (high res version available at www-obs.univ-lyon1.fr/eric.emsellem/preprints/NGC1068_Emsellemetal_final.pdf

Dynamical evolution of AGN host galaxies -— gas in/out-flow rates in seven NUGA galaxies

To examine the role of the host galaxy structure in fueling nuclear activity, we estimated gas flow rates from several kpc down to the inner few 10 pc for seven nearby spiral galaxies, selected from the NUclei of GAlaxies sample. We calculated gravitational torques from near-infrared images and determined gas in/out-flow rates as a function of radius and location within the galactic disks, based on high angular resolution interferometric observations of molecular (CO using Plateau de Bure interferometer) and atomic (H I using the Very Large Array) gas. The results are compared with kinematic evidence for radial gas flows and the dynamical state of the galaxies (via resonances) derived from several different methods. We show that gravitational torques are very efficient at transporting gas from the outer disk all the way into the galaxies centers at ~100 pc; previously assumed dynamical barriers to gas transport, such as the corotation resonance of stellar bars, seem to be overcome by gravitational torque induced gas flows from other nonaxisymmetric structures. The resulting rates of gas mass inflow range from 0.01 to 50 M⊙ yr^(–1) and are larger for the galaxy center than for the outer disk. Our gas flow maps show the action of nested bars within larger bars for three galaxies. Noncircular streaming motions found in the kinematic maps are larger in the center than in the outer disk and appear to correlate only loosely with the in/out-flow rates as a function of radius. We demonstrate that spiral gas disks are very dynamic systems that undergo strong radial evolution on timescales of a few rotation periods (e.g., 5 × 10^8 yrs at a radius of 5 kpc), due to the effectiveness of gravitational torques in redistributing the cold galactic gas

Recent Developments

Context. Tracing nuclear inflows and outflows in active galactic nuclei (AGNs), determining the mass of gas involved in them, and their impact on the host galaxy and nuclear black hole requires 3D imaging studies of both the ionized and molecular gas. Aims. We map the distribution and kinematics of molecular and ionized gas in a sample of active galaxies to quantify the nuclear inflows and outflows. Here, we analyze the nuclear kinematics of NGC 1566 via ALMA observations of the CO J:2-1 emission at 24 pc spatial and ∼2.6 km s−1 spectral resolution, and Gemini-GMOS/IFU observations of ionized gas emission lines and stellar absorption lines at similar spatial resolution, and 123 km s−1 of intrinsic spectral resolution. Methods. The morphology and kinematics of stellar, molecular (CO), and ionized ([N II]) emission lines are compared to the expectations from rotation, outflows, and streaming inflows. Results. While both ionized and molecular gas show rotation signatures, there are significant non-circular motions in the innermost 200 pc and along spiral arms in the central kpc (CO). The nucleus shows a double-peaked CO profile (full width at zero intensity of 200 km s−1), and prominent (∼80 km s−1) blue- and redshifted lobes are found along the minor axis in the inner arcseconds. Perturbations by the large-scale bar can qualitatively explain all features in the observed velocity field. We thus favor the presence of a molecular outflow in the disk with true velocities of ∼180 km s−1 in the nucleus and decelerating to 0 by ∼72 pc. The implied molecular outflow rate is 5.6 M⊙ yr−1, with this gas accumulating in the nuclear 2″ arms. The ionized gas kinematics support an interpretation of a similar but more spherical outflow in the inner 100 pc, with no signs of deceleration. There is some evidence of streaming inflows of ∼50 km s−1 along specific spiral arms, and the estimated molecular mass inflow rate, ∼0.1 M⊙ yr−1, is significantly higher than the SMBH accretion rate (ṁ = 4.8 × 10−5 M⊙ yr−1)

Atomic Hydrogen Properties of AGN Host Galaxies: HI in 16 NUclei of GAlaxies (NUGA) Sources

We present a comprehensive spectroscopic imaging survey of the distribution and kinematics of atomic hydrogen (HI) in 16 nearby spiral galaxies hosting low luminosity AGN, observed with high spectral and spatial resolution (resolution: ~20 arcsec, 5 km/s) using the NRAO Very Large Array (VLA). The sample contains a range of nuclear types, ranging from Seyfert to star-forming nuclei and was originally selected for the NUclei of GAlaxies project (NUGA) - a spectrally and spatially resolved interferometric survey of gas dynamics in nearby galaxies designed to identify the fueling mechanisms of AGN and the relation to host galaxy evolution. Here we investigate the relationship between the HI properties of these galaxies, their environment, their stellar distribution and their AGN type. The large-scale HI morphology of each galaxy is classified as ringed, spiral, or centrally concentrated; comparison of the resulting morphological classification with AGN type reveals that ring structures are significantly more common in LINER than in Seyfert host galaxies, suggesting a time evolution of the AGN activity together with the redistribution of the neutral gas. Dynamically disturbed HI disks are also more prevalent in LINER host galaxies than in Seyfert host galaxies. While several galaxies are surrounded by companions (some with associated HI emission), there is no correlation between the presence of companions and the AGN type (Seyfert/LINER).Comment: 54 pages, 7 figures, accepted for publication in AJ. The full-resolution version is available at http://www.mpia.de/homes/haan/research.htm

A Deep Chandra ACIS Study of NGC 4151. II. The Innermost Emission Line Region and Strong Evidence for Radio Jet-NLR Cloud Collision

We have studied the X-ray emission within the inner 150 pc radius of NGC 4151 by constructing high spatial resolution emission line images of OVII, OVIII, and NeIX. These maps show extended structures that are spatially correlated with the radio outflow and optical [OIII] emission. We find strong evidence for jet--gas cloud interaction, including morphological correspondences with regions of X-ray enhancement, peaks of near-infrared [FeII] emission, and optical clouds. In these regions, moreover, we find evidence of elevated NeIX/OVII ratios; the X-ray emission of these regions also exceeds that expected from nuclear photoionization. Spectral fitting reveals the presence of a collisionally ionized component. The thermal energy of the hot gas suggests that >0.1% of the estimated jet power is deposited into the host interstellar medium through interaction between the radio jet and the dense medium of the circum-nuclear region. We find possible pressure equilibrium between the collisionally ionized hot gas and the photoionized line-emitting cool clouds. We also obtain constraints on the extended iron and silicon fluorescent emission. Both lines are spatially unresolved. The upper limit on the contribution of an extended emission region to the Fe Kalpha emission is <5% of the total, in disagreement with a previous claim that 65% of the Fe Kalpha emission originates in the extended narrow line region.Comment: Accepted for publication in ApJ. 28 pages, 9 figure

A Deep Chandra ACIS Study of NGC 4151. I. the X-ray Morphology of the 3 kpc-diameter Circum-nuclear Region and Relation to the Cold Interstellar Medium

We report on the imaging analysis of 200 ks sub-arcsecond resolution Chandra ACIS-S observations of the nearby Seyfert 1 galaxy NGC 4151. Bright, structured soft X-ray emission is observed to extend from 30 pc to 1.3 kpc in the south-west from the nucleus, much farther than seen in earlier X-ray studies. The terminus of the north-eastern X-ray emission is spatially coincident with a CO gas lane, where the outflow likely encounters dense gas in the host galactic disk. X-ray emission is also detected outside the boundaries of the ionization cone, which indicates that the gas there is not completely shielded from the nuclear continuum, as would be the case for a molecular torus collimating the bicone. In the central r<200 pc region, the subpixel processing of the ACIS data recovers the morphological details on scales of <30~pc (<0.5") first discovered in Chandra HRC images. The X-ray emission is more absorbed towards the boundaries of the ionization cone, as well as perpendicular to the bicone along the direction of a putative torus in NGC 4151. The innermost region where X-ray emission shows the highest hardness ratio, is spatially coincident with the near-infrared resolved H_2 emission and dusty spirals we find in an HST V-H color image. The agreement between the observed H_2 line flux and the value predicted from X-ray-irradiated molecular cloud models supports photo-excitation by X-rays from the active nucleus as the origin of the H_2 line, although contribution from UV fluorescence or collisional excitation cannot be fully ruled out with current data. The discrepancy between the mass of cold molecular gas inferred from recent CO and near-infrared H_2 observations may be explained by the anomalous CO abundance in this X-ray dominated region. The total H_2 mass derived from the X-ray observation agrees with measurement in Storchi-Bergmann et al.Comment: 33 pages, 9 figures and 2 table

A Deep Chandra ACIS Study of NGC 4151. III. the Line Emission and Spectral Analysis of the Ionization Cone

This paper is the third in a series in which we present deep Chandra ACIS-S imaging spectroscopy of the Seyfert 1 galaxy NGC 4151, devoted to study its complex circum-nuclear X-ray emission. Emission features in the soft X-ray spectrum of the bright extended emission (L[0.3-2keV]~10^40 erg/s) at r>130 pc (2") are consistent with the brighter OVII, OVIII, and NeIX lines seen in the Chandra HETGS and XMM-Newton RGS spectra below 2 keV. We construct emission line images of these features and find good morphological correlations with the narrow line region clouds mapped in [OIII]5007A. Self-consistent photoionization models provide good descriptions of the spectra of the large scale emission, as well as resolved structures, supporting the dominant role of nuclear photoionization, although displacement of optical and X-ray features implies a more complex medium. Collisionally ionized emission is estimated to be <12% of the extended emission. Presence of both low and high ionization spectral components and extended emission in the X-ray image perpendicular to the bicone indicates leakage of nuclear ionization, likely filtered through warm absorbers, instead of being blocked by a continuous obscuring torus. The ratios of [OIII]/soft X-ray flux are approximately constant (~15) for the 1.5 kpc radius spanned by these measurements, indicating a relatively constant ionization parameter, consistent with the photoionized outflow of a wind-like density profile. Using spatially resolved features, we estimate that the mass outflow rate in NGC 4151 is ~2Msun/yr at 130 pc and the kinematic power of the ionized outflow is 1.7x10^41 erg/s, approximately 0.3% of the bolometric luminosity of NGC 4151.Comment: 45 pages. 18 figures. Accepted to Ap

Dynamical Evolution of AGN Host Galaxies - Gas In/Out-Flow Rates in 7 NUGA Galaxies

To examine the role of the host galaxy structure in fueling nuclear activity, we estimated gas flow rates from several kpc down to the inner few 10 pc for seven nearby spiral galaxies, selected from the NUGA sample (NUclei of GAlaxies). We calculated gravitational torques from near-IR images and determined gas in/out-flow rates as a function of radius and location within the galactic disks, based on high angular resolution interferometric observations of molecular (CO using PdBI) and atomic (HI using the VLA) gas. The results are compared with kinematic evidence for radial gas flows and the dynamical state of the galaxies (via resonances) derived from several different methods. We show that gravitational torques are very efficient at transporting gas from the outer disk all the way into the galaxies centers at ~100 pc; previously assumed dynamical barriers to gas transport, such as the Corotation Resonance of stellar bars, seem to be overcome by gravitational torque induced gas flows from other non-axisymmmetric structures. The resulting rates of gas mass inflow range from 0.01 to 50 solar masses per year and are larger for the galaxy center than for the outer disk. Our gas flow maps show the action of nested bars within larger bars for 3 galaxies. Non-circular streaming motions found in the kinematic maps are larger in the center than in the outer disk and appear to correlate only loosely with the in/out-flow rates as a function of radius. We demonstrate that spiral gas disks are very dynamic systems that undergo strong radial evolution on timescales of a few rotation periods (e.g. 5 times 10^8 yrs at a radius of 5 kpc), due to the effectiveness of gravitational torques in redistributing the cold galactic gas.Comment: 83 pages, 13 figures, accepted for publication in ApJ. The full-resolution version is available at http://www.mpia.de/homes/haan/research.htm

A REVERSE SHOCK in GRB 160509A

We present the second multi-frequency radio detection of a reverse shock in a γ-ray burst. By combining our extensive radio observations of the Fermi-Large Area Telescope γ-ray burst 160509A at z = 1.17 up to 20 days after the burst with Swift X-ray observations and ground-based optical and near-infrared data, we show that the afterglow emission comprises distinct reverse shock and forward shock contributions: the reverse shock emission dominates in the radio band at ≲10 days, while the forward shock emission dominates in the X-ray, optical, and near-infrared bands. Through multi-wavelength modeling, we determine a circumburst density of , supporting our previous suggestion that a low-density circumburst environment is conducive to the production of long-lasting reverse shock radiation in the radio band. We infer the presence of a large excess X-ray absorption column, N H ≈ 1.5 ×1022 , and a high rest-frame optical extinction, A V ≈ 3.4 mag. We identify a jet break in the X-ray light curve at , and thus derive a jet opening angle of , yielding a beaming-corrected kinetic energy and radiated γ-ray energy of erg and erg (1-104 keV, rest frame), respectively. Consistency arguments connecting the forward shocks and reverse shocks suggest a deceleration time of s ≈ T 90, a Lorentz factor of , and a reverse-shock-to-forward-shock fractional magnetic energy density ratio of . Our study highlights the power of rapid-response radio observations in the study of the properties and dynamics of γ-ray burst ejecta. © 2016. The American Astronomical Society. All rights reserved